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Abstract 



We have investigated the superfluid properties of a ring of weakly interacting and degenerate ID Bose gas at 
thermal equilibrium with a rotating vessel. The conventional definition of superfluidity predicts that the gas has 
a significant superfluid fraction only in the Bose condensed regime. In the opposite regime, it is found that a 
superfluid behaviour can still be identified when the probability distribution of the total momentum of the gas 
has a multi-peaked structure, revealing unambiguously the existence of 'superfluid' supercurrent states that did 
not show up in the conventional definition of superfluidity. 



Cet article etudie, dans le regime d'interaction faible, les proprietes superfluides d'un gaz de Bose unidimensionnel 
confine sur un anneau et a l'equilibre thermodynamique dans un referentiel tournant. La definition habituelle de la 
superfluidite predit que ce gaz a une fraction superfluide appreciable seulement s'il est aussi un condensat de Bose- 
Einstein. Dans le regime non condense, nous trouvons cependant qu'il est possible d'identifier un comportement 
superfluide en considerant la distribution de probabilite de l'impulsion totale du gaz : il existe un regime ou cette 
distribution comporte plusieurs pics bien separes, ce qui demontre l'existence de super-courants superfluides qui 
passent inapercus dans la definition habituelle de la superfluidite. 
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Resume 



1. Introduction 



Is the weakly interacting ID Bose gas superfluid ? 

Diverging answers can be found in the literature. The Landau criterion [1] gives a positive answer, since the 
dispersion of elementary excitations is linear at low momenta. In the thermodynamic limit, it is argued in [2] that, 
on the contrary, the ID Bose gas cannot be superfluid at finite temperature: it does not experience any phase 
transition, and the field correlation function vanishes exponentially at large distances, rather than with a power 
law. A calculation based on the Bethe ansatz, with some additional assumptions on the accessible many-body 
states, concludes that superfluidity is possible [3]. As shown in [4,5] one of the subtleties of the issue is that 
there are actually different definitions of superfluidity, one of them based on a static (that is thermal equilibrium) 
property of the system, the other one involving a dynamical response of the system. In ID, these two definitions 
are found to dramatically differ in the thermodynamic limit [5]. 

In this paper, we restrict to the strict thermal equilibrium regime, in a case where the gas can exchange momentum 
with a rotating vessel with walls that are smooth, at least at the macroscopic scale. Note that this differs from the 
usual stirring procedures used in experiments with condensates, where a macroscopic rotating defect is applied 
[6,7,8,9]. We investigate the superfluid properties of the quantum gas in various limiting cases, from the ideal 
Bose gas to the weakly interacting Bose gas with weak density fluctuations, where the Bogoliubov approximation 
applies. We also consider an exactly solvable classical field model that allows also to study the interacting case with 
large density fluctuations. The applicability domain of this model has some overlap with the one of the Bogoliubov 
approximation, see fig. I, and in this overlap domain the two approaches give coincident results. Investigations 
are performed by considering not only the mean momentum of the rotating gas, but also the whole probability 
distribution of the total momentum: this allows to reach a much deeper physical understanding of the static 
aspects of the problem. 
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Figure 1. Scheme of the different approaches used in this paper and the corresponding applicability domains. The three temperature 
scales on the axis differ by a factor (hp/rag) 1 ^ 2 ^> 1 in the weakly interacting regime, 'fluct.' stands for 'density fluctuations'. 



2. General considerations 



2.1. The physical model 

Consider a one-dimensional Bose gas as described in a second-quantization approach by the Hamiltonian: 

fc2 rL <j2 r 

■«() = -—/ dx^(x) — -*(a;) + | dx¥ {x)¥ {x)il{x)ii{x). (I) 
2m J ox 2 J 



The \&(a;) and ~&^(x) operators are respectively the destruction and creation operators for a boson at point x. 
They obey standard bosonic commutation rules [^/(x),^ (x')] = 5(x — x'). The spatial coordinate x runs on a 
ring of length L with periodic boundary conditions, m is the atomic mass, and the strength of local interactions 
is quantified by g. We shall restrict in this paper to the repulsive and weakly interacting case so that we impose 
< g <C h 2 p/m [10] where p is the mean density. 

The Hamiltonian (1) is a good description of a Bose gas in a cylindrically symmetric toroidal trap [11] provided 
(i) the transverse trapping frequencies in the torus are much larger than both the temperature and the interaction 
energy per particle, hui PtZ >• fcsT, pg, and (ii) the radius L/2-r of the torus is much larger than the width of the 
transverse harmonic oscillator ground state which allows to neglect curvature effects in the kinetic energy. In this 
regime, the system is effectively ID with periodic boundary conditions. 

The gas is assumed to be at thermal equilibrium in an uniformly rotating frame, which mimics the presence of a 
rotating vessel containing the fluid; transfer of angular momentum from the vessel to the fluid is assumed to be 
possible, so that thermal equilibrium between the two can be attained. Rotation is then described by the following 
additional term in the Hamiltonian: 

Wrot = -V TO tP, (2) 

where v rot = QL/2-k is the rotation velocity and P is the total momentum operator of the gas: 

P = -ih J dx¥(x)^!{x). (3) 
Notice that the total momentum operator P gives the total mass current in the laboratory frame [12]. 



2.2. Definition of the normal fraction of the gas 



We define the generalized normal fraction /„ of the gas as: 

fn = AT {P) , (4) 

where the expectation value of the total momentum operator is taken in thermal equilibrium and N is the total 
number of particles. The usual normal fraction of the gas is the limit /„ of /„ for v Iot —> 0, and the corresponding 
supcrfluid fraction is 1 — f° [4]. For a rigid body /„ = 1, which means that it is at rest in the rotating frame: the 
fluid is completely dragged by the walls of the rotating vessel. On the other hand /„ = for a pure superfluid: 
even if the vessel is rotating, the fluid remains at rest in the laboratory frame. 

In the following sections, we shall study in detail the behaviour of f n for a weakly interacting one-dimensional 
Bose gas in different temperature and density regimes. The temperature is assumed to be always much larger 
than the spacing of single particle levels: 

that is the size L of the system is assumed to be much longer than the thermal de Broglie wavelength A: 



/ 27rft 2 

For our finite-size system, only velocity boosts v which are integer multiples of the characteristic velocity: 

are allowed 1 . Galilean invariance under such boosts implies that if ip is an eigenstate of (1) of energy E and 
momentum P, the boosted state ip' is also an eigenstate of energy E' = E + vP + Nmv 2 /2 and momentum 



1 The effect a velocity boost of velocity v on the many-body wavefunction ift(xi, . . . , xjy) is 

*//(xi,...,x N ) =ip(x 1 ,...,x N )e zm,, '£i x i ' h . (8) 
Periodic boundary conditions therefore impose rav — 2irhM / L , M being an integer. 
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P' — P + Nmv. Provided the mean number of particles N is kept constant, the function P(v rot ) giving the mean 
momentum (P) as a function of velocity v TOt therefore satisfies a periodicity condition of the form: 

P(v TO t +vi) = P(urot) + Nmvi. (9) 

Because of the symmetry under spatial inversion, P(v T ot) is an odd function of v ro t- This property, combined with 
(9), implies that P(v\/2) = mNvi/2 and therefore /„(t>i/2) = 1. For this reason, in the following we shall restrict 
the definition (4) of the generalized normal fraction to the velocity range v TOt G [— vi/2, vi/2]. 

2.3. Probability distribution of the total momentum P 

In the previous subsection we have introduced the concept of normal fraction /„ of the gas in terms of the 
expectation value of the total momentum operator P. We shall extend our analysis by considering not only the 
average value of P, but rather the complete probability distribution p(P) , which gives the probability for the total 
momentum to be equal to some given value P. Notice that this probability distribution p(P) is totally different 
and distinct from the usual momentum distribution, which gives instead the mean number of particles in each 
momentum state. 

As usual, the first step for the determination of p(P) is the calculation of the corresponding characteristic function 
5(C): 

5(0 = (e <P ). (10) 
The probability distribution p(P) is then obtained as the Fourier transform of g: 

p(P) = j^e- lCP g(0- (11) 

The slow rotation v ro t — * limit of the normal fraction can be related to the variance of p(P) in the non-rotating 
Vrot = system. As the total momentum P commutes with the Hamiltonian TCo, the density matrix of the fluid 
in the canonical ensemble at a given inverse temperature (3 = 1/fcsT can be expanded for small v TO t as: 

p = e -/3(H -Pv Tot ) _ e -0H o ^ + p Vmt p} . ( 12 ) 

This expression can be used to calculate the expectation value of the momentum P and then the small velocity 
limit /° of the normal fraction. 

/° = lim /„= I ' (13) 
urot— >o m Kb -I jV 

This relation also holds in the grand-canonical ensemble, as P commutes with the number of particles and as the 
chemical potential /i varies only to second order in v TOt for a fixed mean number of particles. It is easy to see 
that for a Boltzmann gas of N distinguishable and non-interacting particles, the equipartition theorem of classical 
statistical mechanics implies that 

(P 2 ) = <P?> = mNk B T, (14) 

i 

where pi is the momentum of the i-th particle. The system is therefore totally normal f° = 1. 



3. Non-interacting gas 



The present section is devoted to a study of the rotational properties of a non-interacting gas. In the grand- 
canonical ensemble, the population of each one-particle state is described by the usual Bose distribution: 

nfe = TP^T) - 1 ' (15) 



where the momentum k is quantized as usual as: 

2-nh 



n, (16) 
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n being an integer, and the single-particle energy e k (v T ot) in the rotating frame at v rot is equal to: 



The mean number of particles is 



and the mean momentum is 



h 2 k 2 

e k {v TOt ) = — hkv Iot . (17) 

2m 



(P)=^hkn k . (19) 

k 

The normal fraction /„ is then immediately obtained from its definition (4). Its zero- velocity value /„ could also 
be obtained from (13) by taking into account the fact that for an ideal gas one has: 

(n\) = 2 {n k } 2 + {n k } . (20) 

As for Vrot = 0a Bose- Einstein condensate can only appear in the k = mode, the prediction for /° is not affected 
by the non-physical grand-canonical fluctuations in the condensate mode. 

3.1. Non degenerate gas 

In the limit f3/j> — > — oo, the occupation of all single particle modes is much smaller than unity and the Bose 
distribution (15) can be approximated by a Maxwell-Boltzmann law of the form: 

n k =e~^"-^. (21) 

As n k is a slowly varying function of k (the condition (5) is always assumed), one can replace the sum over k in 
(18) and (19) by integrals ^ J dk. By switching to the integration variable k' — k + mv/h, it is then immediate 
to see that: 

{P)=mNv Iot , (22) 

which means that the fluid is completely normal and no superfluid fraction is present. Notice that this conclusion 
does not depend on the specific choice of the Maxwell-Boltzmann distribution (21). It rather depends on the facts 
that (i) the occupation number of the mode of wavevector k is a function of hk — mv Iot only, and (ii) sums can be 
replaced by integrals. Notice that the condition (i) no longer holds in the presence of interactions, which allows 
e.g. for superfluid behaviour of the 3D Bose gas even in the thermodynamical limit. For the ideal Bose gas, as we 
shall see in the next subsection, condition (ii) is violated in the Bose-condensed regime. 

3.2. Classical field approximation 

If one assumes that the temperature is larger than the absolute value of the chemical potential fcsT >• \fj,\, which 
corresponds to the limit f3fj, 0~ , the Bose distribution (15) can be approximated by the classical field one: 

n k = . (23) 

£fc - (J, 

Under this approximation, analytical results can be obtained for the normal fraction /„. 
The mean density N can be written as an infinite sum as: 

N = V ^1 - mL2fcaT V - (24) 

where v = v TOt /vi is the rescaled velocity and v 2 — —mL 2 fj,/2ir 2 h 2 . By applying the Poisson summation formula: 

£/(n)=£/(2™), (25) 

ngZ ngZ 
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where f(x) is an arbitrary function and f(k) = J dxf(x)exp(—ikx) its Fourier transform, one is led to the final 
expression: 

_ mL 2 k B T 1 sinh(27rfo) . . 

h 2 2ttu cosh(27rz/ ) — cos(27to) ' 

where z/o is defined as Vq — v 2 — v 2 2 . By applying the same Poisson summation formula (25) to the mean 
momentum (P), one obtains: 

, nl A7 mLksT sin(27ri5) 

(P) = Nmv -2 — i ' — — . (27) 

N ' ft cosh(27ri/ ) - cos(2tto) v y 

Notice the periodicity of iV and (P) as functions of the reduced velocity v. This in agreement with the general 
result (9). From (26) and (27) we immediately obtain the generalized normal fraction: 

r = (P) =1 _!± Sirl ( 27r ") (2g) 

Nmv v sinh(27rz/ ) ' 

The physics of the non-interacting classical field is determined by the dimensionless parameter v and the rotation 
velocity v. For a given rotation velocity v Iot , the crossover from v 3> 1 to v <C 1 corresponds to the Bose 
condensation in the finite system. In the high density v — > limit where the gas is fully Bose-condensed, one has 
vq — > iv and the normal fraction (28) tends to zero. On the other hand, for v — » oo, va — » oo and the normal 
fraction consequently tends to 1. 

In the absence of rotation v = 0, Bose condensation in the finite system occurs for [13]: 

ksT ~ fcsTBEC — ( 29 ) 

that is when the coherence length l c — p\ 2 /2n of the gas becomes larger than the size L of the system. Even in 
this zero-velocity case, the normal fraction 

^ = 1-t§4 (30) 
sinh(27rz/o) 

differs from the non-condensed one / nc : 

, _ N - (nfc=p) _ cosh(27T^ ) - 1 
N ■kvq sinh(27rz/o) 

In particular, in the limit of a well-established Bose- Einstein condensate i/j « 1, one has the peculiar relation: 

fn = 2/nc, (32) 

which states that the normal fraction is twice as large than the non-condensed one, and, therefore, that the 
condensate is not completely superfiuid. This prediction is qualitatively opposite to what is found in liquid 4 He, 
which at low temperatures is totally superfiuid even if the condensate fraction is only of the order of 10% [14,15]. 



3.3. Probability distribution of the total momentum 



For the ideal Bose gas, the total momentum is simply written in terms of the single-particle occupation numbers 
n k as P — ftkrik, so the characteristic function g p is equal to: 

9(o = u- ^ nkn " ) = n {^ hk nk ) = n T^k^w • ( 33 ) 

k k 

which has a factorized form over single particle states. In the classical field approximation, this infinite product 
can be calculated exactly, as shown in section 5. 

The probability distribution p(P) without the classical field approximation is easily obtained by numerical Fourier 
transform of (33) (fig. 2). For a vanishing v TOt , the distribution p(P) is simply narrowed as the temperature is 

2 For v 2 < v 2 , which happens when L < [1 — cos(27rv)]ph 2 /(niksT), vq is a purely imaginary quantity and the analytical 
continuation of (26) has to be taken. 
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Figure 2. Probability distribution p(P) of the total momentum P for an ideal Bose gas of 1000 particles. Rotation velocity Vrot — 
(solid lines), v Tot /vi — 0.125 (dashed lines). Temperatures T j (Nmvf/2) — 0.05, 0.125, 0.25 (from top to bottom: the corresponding 
normal fractions are respectively: f n — 0.3,0.64,0.92). 

decreased and no additional structure appears. An analytical proof of the fact that p(P) has a single maximum 
for «rot = is given in the appendix A in the classical field approximation. Superfluidity effects can therefore be 
seen only as a more rapid shrinking of the distribution as compared to the case of a normal gas (cf. (13)). For 
a slow but finite rotation speed v TOt < vi/2, the flow of particles simply appears as a shift-like distortion in the 
probability distribution 3 . 



4. Interacting gas I: multi-valley Bogoliubov approach 

The present section is devoted to the calculation of the rotational properties of a weakly interacting Bose gas 
at temperatures low enough for density fluctuations to be weak. In this regime, the Bogoliubov theory can be 
used to describe the system. Since we are dealing with one-dimensional systems in which no long-range order is 
present for the phase of the Bose field, the extension of the Bogoliubov theory to quasi-condensates with only a 
finite-range order should be used, as discussed in [16]. For simplicity we shall use the version of the method in the 
canonical ensemble. 

Quantitatively, the weakly interacting gas condition can be written in terms of the density p — N/L and the 
healing length £ as [10]: 




(34) 



while the weakness of density fluctuations Ap/ p <JC 1 requires [16]: 




(35) 



In the weakly interacting regime (34), T d f is always much lower than the temperature 

Tdcg = 2-Kh 2 p 2 /m 



(36) 



for quantum degeneracy. 



3 Calculations are done here for the ideal Bose gas in the grand canonical ensemble. In case where a condensate is present, this 
ensemble will give predictions for p(P) in agreement with the canonical ensemble only if the condensate is in the mode k — 0. 
Otherwise the non-physical grand canonical fluctuations in the condensate mode will dramatically affect p(P). 



4.1. Mean total momentum in the Bogoliubov approximation 



The preliminary step to the Bogoliubov approach is to use the pure quasi-condensate approximation and identify 
the quasi-condensate wavefunction <t>o(x) as a local minimum of the Gross-Pitaevskii energy functional: 



E G pW\ = / dx 



77— \d x (/>\ 2 + ^|0| 4 + itiVrot<t>* d x q 



2m 1 rl 2 

A local minimum has to be a stationary state of the Gross-Pitaevskii equation 



n4> - 



■— d\ + iVff|0| 2 + ihv TO td x 
2m 



(37) 



(38) 



Here we take for stationary states the usual form of plane waves: 



Mx) = -±=e k °*. (39) 

The periodic boundary conditions impose that the momentum is quantized as fco = 2itwo/L, wo being an integer 
called the winding number. The corresponding velocity vq = rao"i will be called in the following quasi-condensate 
velocity. The chemical potential is then 

(i — gp + ^mvo — mvo v ro t- (40) 
The resulting energy of the pure quasi-condensate is then: 

E G p [vrot, vo] = iV^imuo -m« u r otj + ^Npg. (41) 

Under the low-temperature and weak-interaction conditions stated above, the fluctuations of the Bose field around 
the pure quasi-condensate field can be described by the Bogoliubov Hamiltonian [16,17]: 

Hb os = E G p[vrot,v } - Vk + tkblbk, (42) 

where the bk,b\ operators are the usual (bosonic) destruction and creation operators for the Bogoliubov quasi- 
particles. Their dispersion relation efc(«tot,«o) is given by: 

£fc(Vrot, v ) = el + hk(v - v ro t) (43) 

in terms of the usual dispersion relation e° k for a system at rest (vo = t^ot = 0) : 

The Vk coefficients are defined as: 

2m T ^> J y 2m 

The fact that the Vk do not depend on vq — v TOt is a consequence of Galilean invariance. One has also to calculate 
the total momentum operator (3): 

P = Nmv + Y hk b\b k . (46) 

The fact that the total number of particles enters in the first term of the righthand side of this equation (rather 
than the number of particles in the quasi-condensate or the number of particles in the superfluid fraction) is again 
a consequence of Galilean invariance. A derivation of 7-feog and P is given in the appendix B. 

The quasi-condensate mode (j>o{x) is a local minimum of the Gross-Pitaevskii energy functional if all the e k are 
positive, which guarantees the thermodynamical stability of the gas. This condition can be reformulated as: 

r e° 1 

I «o - Vrot I < min fe7 ,o [^|J • (47) 

When the length L is larger than the healing length £, this criterion reduces to the usual Landau criterion 
which states that the flow velocity as measured in the moving frame has to be smaller than the sound velocity 
c = (pg/m) 1 / 2 of the gas at rest. 
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A key assumption of our multi- valley Bogoliubov approach is that the density matrix of the system is a statistical 
mixture of states of different winding numbers w and thus quasi-condensate velocities v — w vi. Each value 
of vo corresponds to a separate minimum - a valley - of the Gross-Pitaevskii energy functional (37) and the 
fluctuations around each of them are treated within the Bogoliubov approximation previously discussed. 

The expectation value of any observable, e.g. the momentum P, is then expressed as an average over the different 
valleys: 

Obviously, only stable states satisfying (47) are to be taken into account. The Z VQ are the contributions of the 
different valleys to the total partition function: 

Z vo = Tr[e-" n U = A e"^(— TT (49) 

where A is a overall factor which does not depend on wo and therefore drops out from all the calculations that 
follow. The total partition function is Zt — v Z V(j . 



The expectation value of P within a given valley (P) vo is obtained as the thermal average of (46) 

exp{— fi[e\ — hk(v TOt — v )]} 



(P) vo = Nmv + £ ,_,ro_iL-,, f ,m- 1 = Nmv ° + Nm ^ ^ ~ ^ < 50 ) 



fe^O 

The function f„ (v) has a very simple physical interpretation: it is the generalized normal fraction of the gas that 
one would predict in a single valley Bogoliubov treatment (i.e. winding number w — 0) for a rotation velocity v. 
Hence the superscript v. 

In the remaining part of this section we shall restrict to the case 

L»£. (51) 
This will provide a considerable simplification of the expression of (P). 

First because this allows to replace in f„(v) the sum over k by an integral. The resulting integral may be calculated 
in the limit of low and high temperatures. If the temperature is sufficiently low: 

k B T«.pg(l- |H| ), (52) 

only the linear part of the spectrum effectively contributes to the integral. This leads to 

_ (k B T) 2 TV 1 / co\ 

JnK) C 3 3mhp(l~vyc 2 ) 2 ' 1 ' 

where c = sj pg/m is the sound velocity. 

For temperatures ksT 3> pg but still smaller than the temperature Tdi at which density fluctuations becomes 
important, a classical field approximation can be performed on the Bose distribution law, which leads to 

Within their own domain of validity, both (53) and (54) predict a single- valley normal fraction much smaller than 
one. 

Let us assume that \v Iot — vo\ <§C c for all the relevant terms in the sum (48). In this case, (50) can be simplified 
as follows: 

(P) vo ~ Nmf^°v Iot + Nmv (l - ft ), (55) 
where f„° is the zero-rotation limit of f„: 

^ = ^°) = ]^^ 2 ^- (56) 

The result (55) can be physically understood in the following terms: for any value of the quasi-condensate velocity 
«o, the normal fraction f„ within the valley is moving at a rotation velocity v TO t, while the superfluid fraction 
1 — fn° moves independently at vo- 
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In this limit \v TOt — vo\ <C c, the expression (49) for the partial partition function for each valley can be simplified 
by expanding the product to second order in v rot — v : 



n i 

11 1 _ e-/3e k (v rotl vo) j-j. -, 

fe^O k=iO 1 



n w 



By inserting this expression into (49), the final expression for the expectation value of the momentum (48) gets 
the simple form: 

P{vrot) = f^°Nmv Iot + N(l - f n °) Y, ff(wd) mvo, (58) 

in which the normal gas always moves at a velocity v TO t and the probability distribution q(vo) for the quasi- 
condensate to have a velocity vo has the Gaussian form: 

q(v ) = N q exp [ - mN ^ k ~f 0) (tW ~ vo) 2 } , (59) 

where N q is the normalization factor. Note that only the fraction 1 — f%'° of atoms which are superfluid in a single 
valley treatment is actually involved in (59). The normal fraction of the gas can be obtained from (58): 

f o _ r ,o + Nm(l R*l g^o ^ [ ~ ¥mN(l - g°) tj] 
fcsT E„ o exp[-i/3miV(l-/^ )^] 

Remarkably this expression coincides with formula (11) of [5] if one takes p s = (1 - K-°)p in [5]. 

We briefly come back to the condition \v Iot — vo\ -C c. For «o = 0, this hypothesis is a direct consequence of (51) 
since v rot can be taken in the interval [— ui/2, vi/2]. What happens if vo S> vi? The maximal accessible value w™ ax 
of «o allowed by the probability distribution (59) corresponds to mvo/2 w k B T/N. u™ ax is much smaller than c 
provided teT< iVpg. As: 

fesT = 

iVp 5 ^ T df L ' 

this is automatically satisfied since we have assumed and T <C Tdf . 



(61) 



Even if the normal fraction f%'° within a single valley is always small f„'° <C 1 in the regime of validity of the 
Bogoliubov approximation, the true normal fraction tends to one as soon as several valleys are thermally populated 

k B T>k B T v = ^Nmv 2 . (62) 

In this case, one can replace the discrete sum over vo by an integral in (58): 

(P)~Nmv Iot , (63) 

finding that the transport of atoms in the presence of a finite v YOt is essentially due to the redistribution of the 
population among the different valleys. In other words, as several different values of the quasi-condensate velocity 
vo are accessible, the presence of a finite rotation velocity causes an imbalance of the relative probabilities q{vo) 
and therefore a net matter flow. 

The Bogoliubov prediction for (P) as function of v TOt is plotted in figure 3 for various values of the temperature. 
In the low temperature regime T T v , only the valley with wo = is generally occupied for « rot < vi/2: 
the full normal fraction is close to the single valley prediction, so that the gas is almost fully superfluid. When 
Vrot approaches «i/2, the valleys wo = and wq = 1 become nearly degenerate: the Gross-Pitaevskii energies 
(41) of a condensate at rest vo = and in the first excited state vo = «i are nearly equal. In this case, even if 
T<T„, both valleys can be populated and their relative weights will be strongly dependent on v TOt , see the step 
function in figure 3. For i; ro t > vi/2, only the valley wo = 1 is occupied. The width of the crossover from a value 
(P) <C Nmvrot to a value Nmv Iot — (P) <§C Nm(vi — v ro t) is of the order of 

In the opposite regime of a temperature on the order of T v , where several valleys are thermally populated, the 
step behaviour of P(v ro t) is smoothed out and the gas become normal, as predicted in (63). 
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Figure 3. Average value of the total momentum P as a function of the rotation velocity v lot for a Bosc gas of 1000 interacting atoms. 
Chemical potential pg/ksTv — 0.1. Temperatures: T jT v — 0.025 (a), 0.1 (b), 0.25 (c). Solid lines: many-valley Bogoliubov predic- 
tions. Dashed lines for (b) and (c): full numerical solution of the classical field model of section 5. For curve (a), fcgT — pg/A <C pg, 
outside the applicability range of the classical field model. 

4.2. Comparison with Quantum Monte Carlo simulation 



The predictions of the many- valley Bogoliubov approach can be verified against a Quantum Monte Carlo calcula- 
tion performed using a recently developed stochastic field technique for the interacting Bose gas [18] . As rotation 
is described by an additional single-particle term —Pv vo t in the Hamiltonian, the numerical simulation does not 
present any further complication with respect to the ones previously performed, e.g. for the study of the statistical 
properties of a Bose- Einstein condensate [19] . 




Figure 4. Average value of the total momentum P as a function of the rotation velocity f ro t for a ID Bosc gas of 42 interacting 
atoms. Chemical potential pg/ksT v — 0.017, T/T v — 0.017. Solid line: many-valley Bogoliubov theory. Points: Quantum Monte 
Carlo calculation. 

In fig. 4 we have compared the result of the Monte Carlo simulation with the prediction of the many-valley 
Bogoliubov theory. As the condition L 3> £ is not fulfilled, a numerical evaluation of (48-50) is necessary. The 
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parameters of the figure correspond to a regime in which only a single valley is thermally populated T <^.T V : the 
agreement is good in the slow- velocity regime v TOt <C V\ as well as close to the point v TOt = Vi/2 where the two 
valleys of winding numbers respectively rao = and wo = 1 become degenerate and the system is in a statistical 
mixture of them. 



4.3. Probability distribution of the total momentum 



If the temperature is sufficiently low for density fluctuations to be small, the many-valley Bogoliubov theory 
introduced in the previous subsection can be used to calculate not only the normal fraction but also the complete 
probability distribution p(P) for the total momentum P. 

Using the expression (46) of the total momentum operator in the Bogoliubov approximation, the characteristic 
function g(Q can be written as: 



_ _L ^2 Z vo e iiNmv ° /e iC ^fc hkn k\ _ _L ^2 e -ii 3Nm ( v o- v 'ot) 2 e i(Nmv a j-j 



Z T ^— ' \ / Z T ^— ' 11 1 _ p-/ 3 [ e fr- fifc Kot-i'o)]piCfifc 

v \ ' v o v „ k 1 e fc e 

(65) 

where n k are now the occupation numbers of the Bogoliubov quasi-particle modes of energies e fe . 

It is apparent in the last factor of (65) that the presence of £ can be reinterpreted as a complex shift of v Iot — v 
by an amount -ifcsT(. If we assume that 

Ifrot - vo - ik B TC\ < c, (66) 
we can simplify (65) in exactly the same way which led to (57): 

~ — ^ e~ pNm{1 ~^' aKvrot ~ vo)2/2 e iCmNvo{1 ~^' 0) e iCmNvlotf "° e - mN fn° k B T < 2 / 2 (g7) 
"0 

The distribution function p(P) for the momentum can then be calculated by Fourier transform of the characteristic 
function (67): 

p (p) _ J_ e -l3rnN(l-f^ a )(v Q -v Iat f/2 e -{P-rnv Q N{l-f^ a )- m v Iot Nf^ a \ 2 /{2 m Nfl- k B T) 

"0 

The structure of (68) is physically transparent and can be summarized as follows: 

- Each valley contributes as a peak of width (fn'° NmkBT) 1 ^ 2 in accord with (13) as applied to a single valley. 

- The peak corresponding to each valley is centered at the value mvoN(l — fn°) + mv T otfn° of the momentum: 
as expected, the (single-valley) superfluid fraction (1 — /^'°) moves at the quasi-condensate velocity vo, while 
the normal one f%'° moves at v TOt . 

- The occupation probability of each valley is proportional to a Gaussian involving the kinetic energy in the 
rotating frame of the (single- valley) superfluid fraction. 

So, for temperatures such that: 

T - « T « ^0 (69) 

Jn 

the probability distribution of momentum p(P) is given by a series of narrow isolated peaks. 

The probability distribution p(P) is plotted in fig. 5 for v Iot = and different values of the temperature. At low 
temperatures T <S T v , only a single valley is populated and p(P) shows a single narrow peak around P = 0. 
In this case, the gas is nearly fully superfluid. For temperatures growing across T/T v w 1, higher valleys start 
to be populated as well, more peaks become visible and the width of each of them gets larger. Correspondingly, 
the superfluid fraction decreases to zero. At sufficiently high temperature, the isolated peaks merge into a broad, 
unstructured, distribution. Note that for all temperatures in the figure, the single-valley normal fraction f%'° 
remains always significantly smaller than unity. 

Conventional superfluidity occurs only when the width of the envelope becomes of the order of the spacing 
between peaks, so that only the one whose quasi-condensate velocity v is closest to the rotation velocity v TOt 
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Figure 5. Multi-valley Bogoliubov prediction for the probability distribution p(P) of the total momentum P for an interacting Bosc 
gas of 1000 particles. Rotation velocity v vot — 0. Chemical potential pg/kgTv — 3.2. Temperature: T/T v — 0.25,0.5,1,2,3 (from 
top to bottom). 
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P/Nrnv, 

Figure 6. Multi-valley Bogoliubov prediction for the probability distribution p(P) of the total momentum P for an interacting 
Bose gas of 1000 particles. Rotation velocity v ro t/vi — (solid), 0.125 (dashed), 0.5 (dotted). Chemical potential pg/ksT v — 3.2. 
Temperature T/T v = 0.5. 

results effectively populated. Otherwise, mass transport occurs in the presence of a finite rotation velocity v Iot as 
a consequence of the redistribution of population between the different valleys more than of the distortion within 
a given valley. This behaviour can be clearly observed in fig.6, where the effect of a finite « rot on the probability 
distribution p(P) is shown: although the shape of each single peak is weakly affected, the envelope function suffers 
a dramatic distortion. 

Let us briefly discuss the validity condition of (66). The validity condition of \v TO t —vo\ «c was already discussed 
in subsection 4.1. What is left is to compare fcsT|£| to c. The narrowest features in p(P) are found to be the peaks 
corresponding to single valley normal fractions, of momentum width (/^'° NrnksT) 1 ^ 2 . As is the Fourier 
transform of p(P) , this results in a maximal value of £ on the order of 
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In the regime UbT pg, f%' is given by the v = limit of (54). We then obtain 

fesrC max 1 



c - (W/ 2 K< 1 (71) 
since we assumed L 3> £. In the regime /cbT <C p<7, we use (53) with v = to obtain 

fc B rc max / fie \ 1/2 

The condition (66) then imposes a lower bound on the temperature, 

k B T > hc/L. (73) 

Since hc/L is the energy of the lowest Bogoliubov excitation in a single valley, this condition is equivalent to a 
non-zero temperature regime in the valley. Note that this condition is automatically satisfied in the multi-valley 
regime T > T v for a weakly interacting Bose gas. 



5. Interacting gas II: classical field model 



The results of the previous section have been obtained on the grounds of a Bogoliubov theory which assumes that 
the density fluctuations are weak; this condition requires the temperature to be sufficiently low T <JC T d f where 
Tdf is given in (35). On the other hand, for temperatures larger than the chemical potential but still lower than 
the degeneracy temperature 

fcsTdcg > k B T > pg, (74) 

a classical field model can be applied. In the present section, we shall discuss the predictions of this approach 
regarding the superfluidity properties of the gas. The actual existence of a temperature range (74) is guaranteed 
by the weak interaction condition (34). In this regime, as one can see in the diagram in fig.l, the applicability 
domains of classical field and Bogoliubov theories have a non- vanishing overlap and, in particular, give coincident 
predictions, as we shall see. 



5.1. The model and its solution 



We generalize to the rotating case the classical field model of [13,20,21]. In this generalization of the model, the 
complex field ip(z) has a grand canonical thermal equilibrium distribution P[ip] proportional to exp(— (3E\4)]) 
where E[ip] is the Gross-Pitaevskii energy functional: 



E[tp] = / dz ^— \d z tp\ 2 + ^|i/>| 4 - ■^«rot^*9 z V - HV>| 2 







(75) 



2m 1 Tl 2 lr ' i 

restricting to the configurations of the complex field obeying the boundary condition ip(0) = ip(L). 

Expectation values of quantum observables are obtained by replacing ip with xp, xji^ with i\>* and then averaging 
over the thermal distribution P[ip]. At this stage, the reader may argue that a classical field thermal distribution 
is expected to lead to divergences in the observables in the absence of an energy cut-off, reminiscent of the 
black-body catastrophe of 19th century. A very fortunate consequence of the ID character of the gas is that the 
classical field model gives finite predictions for the observables relevant for this paper, like the mean density, the 
mean momentum of the gas, the probability distribution of the total momentum of the gas, as we shall see. This 
suppresses the issue of an energy cut-off dependence. 

Calculation of expectation values can be performed exactly in the classical field model: the summation over all 
possible complex paths z — * i>(z) can be viewed formally as a Feynman path integral over trajectories of a single 
quantum particle in 2D, z playing the role of a fictitious time, the real and imaginary parts of ip corresponding to 
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fictitious coordinates x and y. Using in the reverse order the Feynman formulation of quantum mechanics, one can 
map the functional integral over all paths into a Feynman propagator for a fictitious Hamiltonian of a quantum 
particle moving in 2D, here in imaginary time. More details are given in [13,21]. We give here without proof the 
expression of the fictitious Hamiltonian for a rotating system: 

rtj Px+Py, .mV rot 1 . 2 , 2x2 fc a , 1 2w 2 . 2n , na \ 

rot = 2M l ~h~ 2 + V ) - + 2 m ^ot )(x + y ) (76) 

where the fictitious mass is 

Px,p y are the momentum operators of the fictitious particle along x,y and L z — xp y — yp x is the angular 
momentum operator of the particle along z. Note that the fictitious Hamiltonian is not hermitian for v TOt 7^ 0, 
but its anti- hermitian part commutes with its hermitian part, which is indeed rotationally invariant. A numerical 
diagonalization of the hermitian part of H is therefore very simple, as one has to solve a Schrodinger equation for 
the radial part of the eigenfunctions only. The corresponding eigenvectors are labeled by two quantum numbers, 
the angular momentum I £ Z and the radial quantum number n € N, and E n ,i is the corresponding real eigenvalue. 



5.2. Exact expressions for relevant observables 



We give the explicit expression for some useful expectation values. The calculation of the mean density in the 
classical field model is required to determine the chemical potential /j, for a given mean total number of particles. 
The mean density is given by 

{r*P) = (x 2 +y 2 ) q (78) 
where the "quantum expectation value" of any operator O for the fictitious particle is defined as 

Tr \Oe- Ln "rot/ h ] 

<0) « = lt[e-*W»] • (79) 
The characteristic function for the total momentum P — —ih J ip*d z ip is 



M = (£!CP) = Tr[e-^J] (80) 
where we have introduced the complex velocity 

w = Vrot - iC,k B T (81) 

and H w is obtained by replacing v Iot with w in H Vtot . The mean total momentum of the gas is related to the 
derivative of the characteristic function in £ = 0: 

(P) = = Nmvrot - i^^-L{L z ) q (82) 

where we used (78) to obtain the mean total number of particles N. The first order expansion in v rot of this 
expression, when combined with the v Iot — + limit of (4), leads to the exact expression for the standard definition 
of the normal fraction [22] : 

o mk B TL (Lj) q (v rot = 0) 

fi5 • (83) 

Note that this expression makes explicit the fact that one has always /° < 1, which justifies the name of normal 
'fraction'. 

In formula (83) the presence of a non-zero superfluid fraction is related to a non-vanishing expectation value 
{L 2 z ) q . This allows to conclude generally that the superfluid fraction tends to zero in the thermodynamical limit: 
as L — > +00, exp(—LH. Viat /h) becomes proportional to the projector on the ground state of the Hermitian part 
of the fictitious Hamiltonian. As this ground state has a vanishing angular momentum (I — 0), {L 2 z ) q tends 
to zero in the thermodynamical limit. For a non-zero v Iot , one gets similarly that (P) /(Nmv Iot ) — > 1 in the 
thermodynamic limit. One sees that the corresponding critical length scale is L ~ h/SE, where 8E is the energy 
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difference between the first excited state and the ground state of the hermitian part of TC. As the minimal energy 
within a given subspace of angular moment is an increasing function of \l\, 6E is either £ , n=l i=0 — E n=0 j =0 or 
E n=0 ,i=i — E n=0 ,i=o- The lengths corresponding to these two possibilities have been identified in [13], they are 
respectively the correlation length and the coherence length of the bulk. As discussed in [13], the coherence length 
is actually always larger than the correlation length. One then sees very generally that superfluidity in the spirit 
of (4) is exponentially suppressed when the length of the sample greatly exceeds the bulk coherence length. 

The Bogoliubov approach in previous sections of this paper has produced a physical picture in which a superfluid 
behavior can still be identified with valleys even when the standard definition (4) gives a normal fraction close 
to unity. Within the classical field model we can test this prediction in an exact manner, without relying on the 
Bogoliubov approximation. It is useful first to identify the dimensionless parameters on which the classical field 
model actually depends. Let us express the field ip in units of p 1 ^ 2 (where p = N/L is the mean density) and the 
spatial coordinate z in units of 

Lo = -^. (84) 

Note that Lo is on the order of the coherence length of the bulk gas, whatever the value of \ is [13]- O ne then 
realizes that /3E[ip], and therefore the state of the gas, depends only on (i) v la t, the velocity v la t in units of 
ksT I '(ph), (ii) L, the length L in units of Lo, and (iii) on a dimensionless parameter x controlling the interaction 
strength: 

where £ is the healing length such that h 2 /rn£ 2 = pg, T d f is the temperature upper bound (35) required to have 
weak density fluctuations and T v was defined in (62) in the Bogoliubov approach as the temperature lower bound 
to have several valleys populated. 

For a non-rotating gas, « rot = 0, we explore the plane of the two remaining parameters, L and \, by a numerical 
diagonalization of the fictitious Hamiltonians 7i Vrot and Tt w . This gives access to the probability distribution of 
P without approximation, and allows to see in which parameter range this distribution has several peaks. The 
result is plotted in figure (7). As expected, presence of a multi-peaked structure requires a length L larger than 
Lo, otherwise the gas is in the Bose-condensed regime. It also requires a large enough \, that is weak enough 
density fluctuations. The boundary of the multi-peaked domain is studied analytically in the next subsection, 
using a large \ expansion. 

For a rotating gas, we have also performed numerical diagonalizations of Ti. Vlot and TL W , which requires in the case 
of TL W the diagonalization of a non-Hermitian matrix. We have recovered the phenomenon, obtained within the 
Bogoliubov approach, that the peaks of a well resolved multi-peaked structure for p{P) are essentially not shifted 
by the rotation of the vessel, but their amplitudes depend on v Iot . The mean momentum of the gas in the classical 
field model is close to the Bogoliubov prediction in its validity domain, that is for weak density fluctuations, see 
figure (3). 



5.3. Asymptotic expressions 



Analytical results can be obtained in two extreme cases. First, in the ideal Bose gas case, where \ — 0. The 
fictitious Hamiltonians appearing in (80) are then quadratic in the position and momentum operators and can 
be diagonalized exactly. The difference with the usual harmonic oscillator case is that the potential energy term 
MQ, 2 (x 2 + y 2 )/2 is now complex. But one just has to choose for the "oscillation frequency" Q, the determination 
of the square root such that Re 51 > 4 . In this case the usual Gaussian wavefunction oc exp[— mQ(x 2 + y 2 )/2h] 
is a perfectly normalizable "ground state" , and the usual repeated action of the creation operators can be used to 
obtain the "excited states". As a consequence the usual 2D isotropic harmonic oscillator spectrum is recovered, 
£n,i = (2n+ |/| + 1)Q + imwl where n is radial quantum number and I is the angular momentum quantum number. 
The characteristic function can then be calculated exactly as the sum of a geometrical series: 



4 This is not possible if MQ 2 is real negative, which can occur only for £ — and fi + mv ro t 2 /2 > 0. One can then still use the 
formulas to come provided that one uses analytic continuation. 
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Figure 7. From a full numerical solution of the classical field model, domain in the \ ~ L plane where the probability distribution 
of p(P) is multi-peaked. Only peaks higher than 10 -3 times the maximal value of p(P) were considered, and x an d ^ were varied 
in steps of one. 
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cosh(fl Vrat L) — cos(mv Iot L/h) 
cosh(£l w L) — cosh[mL((kBT — iv TO t)/h] 



where the complex oscillation frequencies are such that 
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2m 
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and Q„ 



2m 1 2 



(86) 



(87) 



and w was introduced in (81). The resulting p(P) is shown in the appendix A to have a single maximum, at least 
for «rot = 0. The chemical potential was already calculated in the classical field approximation, using the Poisson 
summation formula, see (26). Using (82) one also recovers the expression (27) for the mean momentum. So, for 
the ideal Bose gas in the classical field model, using the Feynman formula to relate a path integral to the trace of 
an evolution operator is similar to the use of the Poisson summation formula! 

Second, in the large \ limit, where intensity fluctuations of the field become weak, one can obtain asymptotic 
formulas for (/(C)- In this limit, the coherence length Lo is much larger than the healing length £. To simplify the 
calculation, we take v TOt = in a first stage and we restrict to the case of a length L on the order of a few times 
the coherence length and therefore much larger than £ [13]: in this case, the energy differences E n=1 j — i? n=0 ,; are 
much larger than h/L so that, within each subspace of fixed angular momentum I, one can restrict to the ground 
state of the Hermitian part of the fictitious Hamiltonian in the calculation of g(C): 



g(C) oc ^ e - L( - E o,i+ ml ^ k B T )/ h 



lei 



where the normalization factor is obtained from g(0) — 1 and i?o,i is the lowest eigenenergy of the Hermitian 
part of Hw with angular momentum I (remember that i> ro t = here). Using polar coordinates r,9, we write the 
corresponding wavefunction O j as 

AW 



<Po,i(x,y) = fi(r)- 



(2 7 rr) 1 /2 

where r = (x 2 + y 2 ) 1 ^ 2 / p 1 ^ 2 . The purely radial wavefunction then solves the Schrodinger equation 

"^ /i(r) + t/;(r)/;(r)= ^ tM/i(r) 
with an effective potential Ui including a centrifugal term: 

= + \xr A - (A - C 2 /2)r 2 



(89) 



(90) 



(91) 
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where £ = pfK, E ,i = LoE 0t i/h, and where the reduced chemical potential p, — pL^pjiksT) is here close to its 
bulk value calculated for large \ m [13], since L exceeds a few L : 

A = X+ix 1/2 +0(l). (92) 

In the large \ limit, the ground state is deeply localized in the minimum of Ui(r) occurring at a non-zero distance 
r ; from the origin, solution of U'(ri) = 0. One can then expand Ui(r) in a power series around r = rj, include 
the quadratic part in (r — n) 2 in a harmonic oscillator diagonalization, include the cubic, quartic, . . . terms with 
perturbation theory. Treatment of up to quartic terms with second order perturbation theory turns out to be 
sufficient here: 



= U ^ + r i + -fe^ - + (93) 

where u>i = [U^ (ri)] 1 ^ 2 is the oscillation frequency in the harmonic approximation to Ui(r). The position r; can 
be calculated exactly, since rf is found to be the root of a cubic equation; when (93) is plugged into (88), one 
gets a very good approximation for and, by numerical Fourier transform, a very good approximation for the 
probability distribution of P, as compared to the full numerical solution, for x> 1, see figure 8. A more tractable 
expression can be obtained by realizing that typical values of Q 2 and I 2 are of order y}l 2 5 , and by performing a 
systematic expansion of the cubic equation for n and of (93) in powers of x~^^'- 

E , t - IC = const + \{z - l) 2 + - ±(l 2 - C 2 ) 2 + 0( X - 1/2 ) (94) 

where const depends on \ only, not on I or £. The second term in the righthand side of the equation is 0(x 1 ^ 2 ), 
the third and fourth terms are a priori of the same order O(l). In practice the fourth term is typically 0(x~ 1 '" 2 ) 
so we neglect it 6 . We then obtain a Gaussian expressions for g(Q and p(P). 

The previous calculations are immediately generalized to the case of a non-zero velocity: as L > Lq, p + mv Iot 2 /2 
depends weakly on v rot and can be replaced by the bulk value for v TO t = 0; one then has to replace £ 2 by ( 2 — 2iv TOt ( 
in the above calculation of Eq,i, where v TO t, being independent on x, is O(l). We then obtain for <?(C): 



9(C) = A ^ e -M'-C+« r ot) 2 /2 e - £(2 /( 2 * 1/2 ) (95) 



where A is a constant factor such that g(0) — 1. Performing the Fourier transform of this expression and using 
the Poisson summation formula leads to 

p(P) = Be- (p - L ^ )2/i2L) J2 e- ip - 2 ^ 2xl/2/(2i) (96) 

where P = P/(ph). The constant factor B is such that the integral of p(P) over P is equal to unity. Formula (96) 
is very suggestive, as it is simply a Gaussian envelope centered in P = Lv Iot on top of a periodic train of Gaussian 
peaks separated by 2tv in P space. Therefore the function p(P) is multi-peaked if the envelope has a width larger 
than 2-7T and if each Gaussian of the train has a width less than 27r: 

L 1/2 > 2tt and {x U2 /L) 1/2 > 2tt. (97) 

One recovers the conditions (69) using the fact that f° = x~ X/ ^ hi the classical field regime feT > pg, see (54). 

The asymptotically exact expression (96) is successfully compared with the full numerical solution in figure 8. 
One can also compare it to the Bogoliubov prediction (68): the two predictions are found to be identical up to 
higher order terms (see appendix C). As a consequence, all the physical discussion following (68) also holds for 
the classical field model! From (95) one can also calculate the normal fraction f° for v TOt — by taking the second 



5 This can be obtained by trial and error, but also by the fact that, in the Bogoliubov theory, the narrowest structure in the 
probability distribution of P scales as (Z^) 1 ^ 2 , resulting in a maximal value of £ scaling as ■ In the temperature regime 
fcgT > [i, which is the one of the classical field model, f%'° scales as 1/x 1 ^ 2 , sec (54). 

6 After multiplication of (94) by —L/Lq and exponentiation to get g(C), one realises that g(C), considered as a function of is 
approximately a superposition of narrow peaks centered in integer values / and of width oc (Lo/L) 1 ^ 2 , with an envelope of width 
X 1/4 (L /L) 1/2 . As a consequence, for a given /, (-1 is 0((L Q /L) 1/2 ), and the fourth term in the RHS of (94) is 0(/ 2 /x) = 0(x~ 1/2 ). 
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order derivative of (95) with respect to £. Equivalently one may calculate the variance of P from (96) and use 
(13). One gets two equivalent forms: 



f° 

J n 



| L E,9 2 exp(-L eff g 2 /2) 
Lo V exp(-L cff <j 2 /2) 



E qez (2^) 2 e 



2 „-(27rg) 2 /(2i e ff) 



,1/2 



i(i + x- 1/2 ) 2 £ q 



-(2^ 9 )V(2ieff) 



(98) 



where the reduced effective length is L c £f = L(l + x _1 ' 2 )/Lo- The first form immediately shows that f„ tends to 
one exponentially for L 3> Lo. The second form recovers the formula (11) of [5] in a calculation up to first order 
in x 1//2 where one replaces ps of [5] by p(l — X^ 1 ^ 2 )- 

1.5 




Figure 8. Probability distribution of the total momentum P in the classical field model, for x — 200, v ro t — and L j Lo — 20. 
Symbols: full numerical solution. Solid line: analytical expression (96). Dashed line: semi-analytical prediction resulting from the 
use of Ti as given by the exact solution of U'(ri) — (see text). The dashed line and the symbols are almost indistinguishable. 



6. Conclusions 



We have investigated the superfiuid properties of a ring of degenerate and weakly interacting ID Bose gas at 
thermal equilibrium with a rotating vessel at velocity v ro t- Provided the transverse trapping is strong enough, our 
model is a good description of a Bose gas confined in a toroidal trap [11]. 

Using the conventional definition of the superfiuid fraction, which relies on the variance of the total momentum of 
the gas in the limit v Iot — > 0, we find that the gas has a significant superfiuid fraction only in the Bose condensed 
regime, that is when the length of the ring does not exceed the coherence length of the bulk gas l c ~ p\ 2 , p being 
the mean density and A the thermal de Broglie wavelength. 

To investigate more carefully the regime where the length of the ring exceeds the coherence length, we have 
considered the full probability distribution of the total momentum P. We have identified a regime where several 
peaks appear in this probability distribution, each peak corresponding to a quasicondensate in a plane wave state 
with a given winding number, the analog of supercurrents in superconductors. Each supercurrent state exhibits 
some superfiuid behaviour: in presence of a non-zero v ro t, the peaks in the probability distribution of P are indeed 
not shifted. This allows to define a local normal fraction for an individual supercurrent. Quantitatively, we have 
found that the probability distribution of P shows several isolated peaks provided that the length of the ring in 
units of the coherence length does not exceed the inverse of the local normal fraction. 

In this non-Bose-condensed regime, it is obvious that the conventional criterion for superfluidity based on the 
variance of P is sensitive to the envelope of the distribution but does not catch its multi-peaked structure. To 
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get it, a direct measurement of the total current P is required, which, e.g., could be performed by means of 
the technique proposed in [23]: in a slow-light regime, the dielectric susceptibility of the atoms depends on the 
local value of the matter current so that the phase accumulated by light after a round-trip around the ring is 
proportional to the total current P. 



Appendix A. Absence of multiple peaks in p(P) for the ideal gas 



The characteristic function for the ideal gas in the classical field regime is: 

.<o-n Wa *V-U - (AI) 

In the v Iot = case, we can regroup the pairs ±fc and rewrite (A.l) as: 

<0 -S ,M -S ^w ' <A2) 

The Fourier transform Pk{P) of each term g k {() is then: 

Pk(P) = ^ e-"*l"l, (A.3) 

with rjk = (3(e k — fj,)/hk. In particular, Pk{P) has the property of being an even function that is decreasing for 
P > (let's call this property V). As the characteristic function g(£) is the product of the g k (C), the distribution 
function p(P) is the convolution of the p k (P). As the convolution of two functions with the property V gives 
again a function with the property V (a sketch of the proof is given below), we can conclude that the probability 
distribution p(P) for the ideal gas has a single maximum, which is at P — 0. The possibility of multi-peaked 
structures is therefore ruled out for the ideal gas in the classical field regime. 

We can prove that the property V is preserved by convolution operations in the following way. Let pi(P) and 
P2(P) be two arbitrary functions sharing property V. We have to prove that: 

p c (P) = JdP' Pl (P')p 2 (P-P') (A.4) 

i) is an even function. 

ii) is a monotonically decreasing function for P > 0. Let's compute its derivative for P > 0: 

p' c (P) = JdP' Pl (P')p' 2 (P-P') = JdP' Pl (P-P')p' 2 (P') = 

J dP'p' 2 {P') [ P1 (P - P') -pi(P + P')] . (A.5) 

As \P — P'\ < \P + P'\ and pi is a decreasing function of the absolute value of its argument, the integrand is 
negative. This guarantees that p' c (P) < for all P > 0. 



Appendix B. Derivation of the Bogoliubov Hamiltonian and momentum operator 



The Hamiltonian (42) can be obtained either by directly solving the Bogoliubov-de Gennes equations for a moving 
system, or, better, by applying Galilean invariance arguments to the well-known case of a system at rest. The 
eigenstates of a weakly interacting Bose gas at rest (vo = 0) are labeled within Bogoliubov theory by the occupation 
number of the bosonic quasiparticle modes {n k }. Thanks to translational invariance, these eigenstates are also 
eigenstates of the momentum: each quasiparticle carrying a momentum hk, the total momentum of the gas is 
given by: 

P[{n k }] = ^hkn k . (B.l) 

fc#0 



20 



Omitting for the moment the rotation energy —Pv Iot , the total energy is given by: 

E Vtot=0 [{n k }] = E G p[v I0 t = 0, v = 0] - ^ e° k V fc 2 + ^ e£ n fc . (B.2) 

The properties of a moving quasi-condensate at vo 7^ can be obtained from the ones of a quasi-condensate at 
rest by transforming the energy and the momentum via a Galilean transformation of velocity v = vo- As discussed 
in sec. 2.1, the total momentum P' in the moving frame is in fact given by: 

P'[{n k }\ = P[{n k }} + Nmvo = J2 hkn k + Nmv . (B.3) 

Inserting back the rotation term ~P'v rot , the energy E' turns out to be: 

/ 1 2 

E [{n k }} = E Vrot=0 [{n k }] + -Nrnv + (v - v TOt )P[{n k }] = 

= E G p[Vrot,V ] - ^ e k V k + + hk ( v ~ V ia t)] Tl k . (B.4) 

The eigenstates and eigenenergies obtained in this way exactly correspond to the ones of the Bogoliubov Hamil- 
tonian (42). For each state, the total momentum (B.3) agrees with (46). 



Appendix C. Comparison of Bogoliubov and classical field theory for p(P) 



In the regime of weak density fluctuations and a temperature ksT 3> pg we compare the expressions for the 
total momentum probability distribution p(P) obtained by Bogoliubov theory, (68) on one side, and by a large x 
expansion of the classical field model, (96) on the other side. At first glance, (96) looks much simpler and therefore 
different than (68). However one has the identity 

L-\P LS rot ) 2 + L-\P 2.g) Y /2 = L-\l + [P Po)" + L" 1 ^-^ (C.l) 



where 

Lv Iot + 2Tvq X 1/2 



P o = 12 :T ■ (C2) 

Expanding Po up to first order in x 1 / 2 and using the fact that f%'° = X _1 ^ 2 m t ne Bogoliubov theory, one 
finds that P — Po coincides with the expression in square brackets in (68). The factor in front of this expression, 
proportional to l//n'°, recovers Z _1 (l + X^ 2 ) within leading order in x 1 ^ 2 ■ The last term in (C.l) coincides with 
the argument of the first exponential factor in (68) when expanded up to first order in x^ 1 ^ 2 ■ 
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